Depinning of domain walls with an internal degree of freedom 
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Taking into account the coupling between the position of the wall and an internal degree of free- 
dom, namely its phase <p, we examine, in the rigid wall approximation, the dynamics of a magnetic 
domain wall subject to a weak pinning potential. We determine the corresponding force-velocity 
characteristics, which display several unusual features when compared to standard depinning laws. 
At zero temperature, there exists a bistable regime for low forces, with a logarithmic behavior close 
to the transition. For weak pinning, there occurs a succession of bistable transitions corresponding 
to different topological modes of the phase evolution. At finite temperature, the force-velocity char- 
acteristics become non-monotonous. We compare our results to recent experiments on permalloy 
nanowires. 
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Many physical systems comprise different phases which 
coexist and are separated by an interface. Examples 
range from magnetic [1, 2, 3, 4] or ferroelectric [5, 6] 
domain walls (DWs), to growth surfaces [7, 8] or contact 
lines [9] . Common to this large variety of phenomena is a 
macroscopic description within which the interface prop- 
erties are well described by a competition between the 
elasticity, which tends to minimize the interface length, 
and the local potential, whose valleys and hills deform 
the interface so as to minimize its total energy. 

Such interfaces are described by the theory of disor- 
dered elastic systems [10, 11], which explains well their 
static (e.g., roughness at equilibrium, correlation func- 
tions) as well as dynamical features (transient regime, 
response to a field). The existence of a threshold force / c 
below which the system is pinned is a crucial feature of 
the zero-temperature motion of such an interface. When 
/ > / c the velocity v ~ (/ — / c ) is characterized by a de- 
pinning exponent /3, while at finite temperature v ~ 
at / = f c defines the thermal exponent ip. Some pre- 
dicted exponents are in very good agreement with mea- 
surements, e.g., in magnetic [1] or ferroelectric [6] films, 
while discrepancies remain for contact lines [9] or for mag- 
netic wires [12], and one can ask what are the missing in- 
gredients in the description. In particular, it is usual in 
the macroscopic description to specify only the position 
of the interface, discarding, a priori, as irrelevant internal 
structures. Here we investigate how this position couples 
to an internal degrees of freedom, and how this coupling 
is manifested in experiment. 

In magnetic systems, the DW position is gencrically 
coupled to an internal degree of freedom (a spin phase 
ip [35]). An interesting case is the motion of a 180° DW 
in a narrow ferromagnetic thin film, which has been the 
subject of intense experimental study because of its im- 
portance for spintronics [12, 13, 14, 15]. It is known [16], 
in the absence of pinning, that i>(/) increases up to a 
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FIG. 1: (color online) (a) Finite temperature, zero current 
(j B = 0) force- velocity curves (blue, solid) and emf ^(<p) (red, 
dashed) for a DW with pinning. Parameters correspond to 
experiments shown in (b). Adapted from Ref. [13], the data 
(squares) are for small but opposite currents. In the absence 
of pinning (Walker model, dashed blue line) there is no second 
rapid increase in v manifest in the data. 



characteristic 'Walker' force / w above which the veloc- 
ity actually decreases up to values / 3> / w (Fig.l). In 
contrast standard interface theory [10] takes pinning into 
account but not the phase, yielding a monotonous v(f). 
An approach combining these two ingredients is missing. 

In this Letter, we develop such an approach and show 
there are dramatic changes as compared with both the 
Walker picture and standard interface theory. Specifi- 
cally, the DW is pinned up to a force /* above which 
the depinning law is bistable and logarithmic. Even more 
strikingly, as / is increased further, the velocity falls back 
to zero until a second depinning transition occurs (Fig.l). 
This is followed by a cascade of such transitions until fi- 
nally v(f) becomes monotonous. Upon adding the effects 
of a finite temperature, this offers a natural explanation 
of the two peaks of v(f) observed in experiments [13]. 



We consider an uniaxially anisotropic ferromagnetic 
medium with position-dependent magnetization of direc- 
tion fi = (sin#cos<;9,sin6>sin</j,cos£?) with easy 2-axis 
and hard y-axis with respective anisotropy constants K 
and K±. With spin stiffness J, the energy E[Si\ is [17] 

E = \j d3x { J ( v6 2 + sin2 °(^ 1 p) 2 

+ K sin 2 6 + K± sin 2 9 cos 2 <p j 
The corresponding Landau-Lifshitz dynamics reads 



(d t + v s W)n = Slx(H- 



fix (ad t + l3 s v s W)n (1) 



Here, H = —5E/5Q, + / ext and / ext is the external field 
and a is the Gilbert damping, which accounts for dissi- 
pation. The velocity v s [36] is proportional to the spin- 
polarized current density j s and [3 S is the current-induced 
relaxation. The white noise r\ accounts for thermal fluc- 
tuations for temperature T [37] . Below the Walker field 
/w = -^ctK±, and for T = 0, there exists a solution to (1) 
for constant / oxt [16]: 6(x,t) — 2 arctanexp [(x/X — 
Vwt)/t;] with £ = [1 + (Kj_/K) sin tp] ~ x / 2 and constant 
tp(x,t) = \ arcsin / oxt // w . This represents a Neel DW of 
width A = \/ J/K and velocity v w = £/ ox t// w [38]. 

In more general situations, this domain wall solution 
with v w t replaced by the actual position r(t) and with 
this and tp(t) considered as parameters is used as an 
ansatz. In the rigid wall approximation (constant £), one 
obtains the effective equations [16, 18, 19, 20, 21, 22] 



ad t r - d t tp - f3 s v s = fext(r) + m 



ad t ip + d t r - v s 



\K± sin 2tp 



(2) 
(3) 



(A = 1 by choice of length units [39].) We split the ex- 
ternal field f C xt(r) — f — V'(r) into a constant 'depin- 
ning' (or 'tilt') force / and a 'pinning' force —V'(r) de- 
riving from a potential V. The effective thermal noise is 
now [22] (r)i(t)r}j(t')) = 2(hN)- 1 ak B T5(t' - t) 5 l0 where 
N = 2 A A/ a 3 is the number of spins in the DW, of section 
A. For constant fields / > / w , this ansatz reproduces 
very accurately [16] the numerical solution of the bulk 
equation (1). We checked that this result also extends to 
the case of a non-uniform potential landscape [23] . 

Having simplified (1) to (2-3), we further restrict our 
study to the case j s = 0: current effects will be consid- 
ered elsewhere [23]. The potential V(r) should reflect 
the pinning effects of impurities or local variations in the 
couplings. A proper treatment of a realistic disordered 
V(r) is a delicate task and in order to gain insights into 
the full problem we take a periodic V(r) = £ sin Kr [40]. 

We first consider the zero-temperature motion. Before 
embarking into a thorough analysis of (2-3), let us gain 
some insights from simple considerations. At / = 0, the 
wall is pinned in one of the minima of V(r) = - sin nr. 
There exists a characteristic force f c beyond which local 
minima of the tilted potential - sin nr — fr disappear 




FIG. 2: (color online) Phase space trajectories (r, ip) for / 
in the bistable regime (/* < / < / c ). S, U, Hi : 2 are the 
fixed points. Blue and turquoise trajectories converge to the 
attracting limit cycle (yellow). Those in green end at the 
stable fixed point S. Separatrices (red) mark the boundaries 
between the corresponding regions. The repulsive limit cycle 
(dashed yellow) is also a separatrix. 

(here f c = 1). Ignoring the variable <p, the wall would 
start to move at / = f c and acquire a finite mean velocity 
at long times for / > f c , because of damping. But since r 
is coupled to ip, the wall may store enough kinetic energy 
in ip to cross barriers for forces less than / c , hence shifting 
the depinning transition to some /* < / c . For / between 
these values, the system is bistable: depending on the 
initial condition, the wall is either pinned in a minimum 
or slides down in the tilted landscape while tp oscillates 
around its own minimum. Moreover, the periodicity of 
ip can induce an unexpected effect: increasing / makes ip 
cross its own barrier and fall into its next minimum, but 
this has a cost: dissipation increases and p> cannot give 
back enough kinetic energy to r. This intuitive picture 
explains the valley appearing in v(f) (Fig.l) until the 
depinning force injects enough energy to reach another 
regime where both ip and r increase in time. 

The analysis of (2-3) can be put on a firm basis by 
considering the phase space of (r, ip) , which is a torus of 
period 27r/«; in r and 7r in tp. We determine the nature 
of the possible trajectories. These cannot cross in phase 
space, but can meet at the fixed points (steady solutions 
of (2-3)). Trajectories approaching a stable fixed point 
have zero mean velocity v = (dtr), unlike those moving 
along a limit cycle. In the latter case the particle covers 
one spatial period 2t:/k over a period of time t, so that 
the average velocity is given by v — 2it/(kt). We will 
thus determine v(f) through r. 

We remark that (2-3) has no fixed points for / > f c , 
which means that the DW moves with non-zero velocity. 
For / < /c, there are four fixed points of coordinates 
(Fig.2): H x = (r ,0) and H 2 = (-r ,7r/2), which are 
hyperbolic (i.e., with one unstable and one stable direc- 
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tion), S — (— r ,0), totally stable, and U = (r ,7r/2), 
totally unstable; here ro = arccos/ > 0. The decomposi- 
tion of the phase space into different dynamical regimes 
depends on the value of K± , which determines how much 
ip can depart from (see [23] for full details): 

(i) Case of high K±: all trajectories end at S, i.e., 
there are no limit cycles (v — 0). This regime persists 
until / = f c , at which point the pairs (Hi,S), (H 2 ,U) 
merge and give rise to a saddle-node bifurcation [24] . For 
small Sf = f — f c > 0, we have adtr ~ 6f + \r 2 in the 
region r ~ where the trajectory spends most of its time, 
whence r(t) = y/2Sf tanjty Sf/2] and one recovers the 
standard depinning law with exponent j3 = J (Fig.3.b). 

(ii) Case of intermediate K±: there is in addition a 
second critical force /* < f c (Fig.3.b). For / < /* all tra- 
jectories end at S. For /* < / < f c one has bistability: 
depending on the initial condition (Fig.2), trajectories ei- 
ther end at S (v = branch in or move along an 
attracting limit cycle (v > branch). In that case, the 
bifurcation is homoclinic [24] and /* is the value of the 
force for which the trajectory starting in an unstable di- 
rection of Hi returns exactly to Hi (see Fig. 3. a). The 
value of /* depends on global aspects, in contrast to the 
case (i). The depinning law of the v > branch is eval- 
uated as follows (see [23]): the trajectory spends most of 
its time close to the hyperbolic point Hi, of positive Lya- 
punov exponent A so that r(t) ~ (/ — /* )e At . The period 
r thus verifies 2tt/k ~ (/ — f*)e Ar and up to factors: 

vtxIlogC/-/*)!" 1 (4) 

(Hi) Case of smaller K±: here, the particularly novel 
features appear, with a depinning law v(f) characterized 
by a succession of bistable regimes (Fig.4.b) separated 
by regions of zero velocity. An interesting mechanism 
emerges: in general, the first bistable regime is charac- 
terized by cycling trajectories with r advancing and ip 
oscillating within a bounded interval, but with increas- 
ing force, ip will eventually rotate by a whole period of tt, 
and fall into S. At this point there is a collision between 
the stable and unstable limit cycles (of Fig.2), and the 
original type of limit cycle disappears for larger forces re- 
sulting in an intermediate v = valley. Increasing / even 
more, the phase space reorganizes until there appear tra- 
jectories with both r and ip increasing with each period 
of the limit cycle (Fig.4.b). Each bistable regime is gov- 
erned by the same bifurcation as in case (ii), but is now 
also characterized by the number of windings of r and ip 
during r. This striking phenomenon, a topological transi- 
tion, arises from the periodicity of the phase and cannot 
appear for instance in situation where r couples to an 
unbounded variable (e.g., the momentum of a massive 
particle in a periodic potential). But topological transi- 
tions can potentially be found in other systems, e.g., for 
viscously coupled particles in a periodic potential [25], 
described by equations similar to (2-3), although v(f) is 
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FIG. 3: (color online) (a) Sketch of the phase space for the 
coordinates (f, (p), and homoclinic bifurcation. For f < f* 
(red, dashed), trajectories all end at S. For f > f* (black, 
solid), the line starting from the unstable direction of Hi con- 
verges to the limit cycle. At T > 0, the random evolution has 
Gaussian probability distribution of variance oc VT. (b) The 
two types of depinning laws at T = 0. For large K± (green, 
dotted) depinning occurs with = |. For intermediate K± 
(blue, solid) there is a bistable regime of force (/* < / < / c ) 
with a zero- velocity branch and a non-zero one in 1/ log. 

found to be monotonous for the conditions of [25] . 

We finally address the finite temperature dynamics in 
the regime /* < / < f c , of particular interest since 
thermal fluctuations cause the system to forget its ini- 
tial condition, and thus destroy the bistability. Taking 
normal coordinates close to Hi, (Fig. 3) the evolution fol- 
lows dtf = er + f 2 + rj, dtp = —<p (with e oc f c — / > 0). 
Starting from Hi towards the limit cycle and evolving 
with the noisy dynamics, the trajectory comes back to Hi 
with a Gaussian distribution of width oc yT at distance 
d(e) oc e — e* from the separatrix (Fig. 3. a). The mean 
escape time is determined by a competition between the 
large Arrhenius time to escape from the local potential 
trap V(r) — —er 2 /2 — f 3 /3 and the small probability 
~ exp(-d(e) 2 /T) of falling into it [23]: 

Escape ~ exp I — - A J (5) 

Thus for T > 0, the bistability curve is transformed in 
the following manner: the curves v(f, T) all cross at some 
new characteristic force /** (where the polynomial in e 
in (5) has a zero) [41]. For / < /**, the depinning is 
dominated by the escape from the trap while for / > 
/**, v(f, T) approaches the positive branch of the T = 
law (Fig. 4. a). In the limit T — » + , v(f) is monostable 
and discontinuous in /**, in contrast with the T = case. 
Note that Vollmer and Risken [26, 27] have studied the 
dynamics of a massive particle in a periodic potential, for 
T > 0, with results related to those of Fig. 4. a, but with 
an approach limited to that particular problem and only 
valid in the small a regime. In contrast, our approach 
not only displays a non-monotonous v(f) but also allows 



4 



V 


(a) 


T = 






> 




----T -> 0+ 











(b) 


W(r) 


= 1 


W(r) = 


1 




= 


V%>) = 


1 



FIG. 4: (color online) (a) For any finite T there exists a criti- 
cal field /** such that for / < /** the albeit small probability 
of falling into the trap S exceeds that of escape and the veloc- 
ity is small while the inverse is true for / > /**. The curves 
approach the vertical line as T — > and all curves cross at 
/ = /** to first order in T. (b) For smaller K ± (and T = 0), 
successive bistable transitions occur, indexed by the winding 
numbers W(r) and W(ip) of r and <p on one period. 



for a general discussion of what happens in the vicinity of 
a homoclinic bifurcation for any a, in the perspective of 
stochastic differential equations, and is more in the spirit 
of Freidlin and Wentzell [28, 29]. 

Despite the oversimplified features of our model, it of- 
fers a possible resolution of an experimental mystery [13]. 
In the absence of pinning (V = 0), and corresponding to 
the simple Walker breakdown picture it is possible to re- 
produce the first peak of v(f) for parameters similar to 
those of [13] (a = 1.32 x 1CT 2 , K± = 1500 Oe, corre- 
sponding to f w = 9.9 Oe, v w — 200 m/s), but v(f) has 
only a peak without the valley seen in the experimental 
data (Fig.l.b). In contrast, due to the existence of a topo- 
logical transition, for finite V the velocity v falls towards 
zero following a first peak and only rises again for larger 
values the force /. Indeed, simulations [42] of (2-3) repro- 
duce the 'valley' observed in v(f) for similar parameters 
(a = 5 x 10~ 2 , K ± = 1200 Oe, V = 50 Oe, k = 5.15/A). 
We predict in particular the following watermark for the 
topological transition: the second peak should coincide 
with the appearance of non-zero (ip) , measurable through 
the emf [30] | (if) (Fig.l.a) while in the Walker picture 
(p) > immediately for / > f w . 

To summarize, we have shown how the coupling be- 
tween the phase and the position of a rigid wall in ID 
dramatically affects the depinning law, which displays 
bistabilities and an unusual scaling v cx 1/| log(/ — /*)|. 
Due to the periodicity of the phase, there are conditions 
for which different bistable regimes follow one another 
with increasing /, yielding for T > a non-monotonous 
v(f), which might well explain features of recent mea- 
surements [13]. It would be valuable to consider the 
current-driven case (v s ^ 0), of interest in the context 
of spintronics. Moreover, the solitonic ansatz can also 
describe interfaces with non-zero dimension, where the 
interplay between the phase and the elastic deformations 
potentially affects the creep motion and the depinning. 
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